Statistical Properties of Eigenstates in three-dimensional Mesoscopic Systems with 

Off-diagonal or Diagonal Disorder 
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The statistics of eigenfunction amplitudes are studied in mesoscopic disordered electron systems of 
finite size. The exact eigenspectrum and eigenstates are obtained by solving numerically Anderson 
Hamiltonian on a three-dimensional lattice for different strengths of disorder introduced either in the 
potential on-site energy ("diagonal") or in the hopping integral ("off-diagonal"). The samples are 
characterized by the exact zero-temperature conductance computed using real-space Green function 
technique and related Landauer-type formula. The comparison of eigenstate statistics in two models 
of disorder shows sample-specific details which are not fully taken into account by the conductance, 
shape of the sample and dimensionality. The wave function amplitude distributions for the states 
belonging to different transport regimes within the same model are contrasted with each other as 
well as with universal predictions of random matrix theory valid in the infinite conductance limit. 
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The disorder induced localization-delocalization (LD) 
transition in solids has been one of the most vigorously 
pursued problems in condensed matter physics since the 
seminal work of Anderson. 1 In thermodynamic limit, 
strong enough disorder generates a zero-temperature crit- 
ical point in d > 2 dimensions 2 as a result of quan- 
tum interference effects (in d < 2 even weakly disor- 
dered metal turns into an Anderson insulator for suffi- 
ciently large sample size). Thus, research in the "pre- 
mesoscopic" era 3 was mostly directed toward the view- 
point provided by the theory of critical phenomena. 4 The 
advent of mesoscopic quantum physics 5 has unearthed 
large fluctuations, induced by quantum coherence and 
randomness of disorder, 6 of various physical quantities 7 
(e.g., conductance, local density of states, current relax- 
ation times, etc.), even well into the delocalized phase. 
Thus, complete understanding of the LD transition re- 
quires to examine full distribution functions of relevant 
quantities. 8 Especially interesting are the deviations of 
their asymptotic tails, a putative signature of incipient 
localization, 7 from the (usually) Gaussian distributions 
expected in the limit of infinite dimensionless conduc- 
tance g = G/Gq (Gq = 2e 2 /h is the conductance quan- 
tum). 

This paper presents the study of such type — numerical 
computation of the statistics of eigenfunction amplitudes 
in finite-size three-dimensional (3D) nanoscale (com- 
posed of ~ 1000 atoms) mesoscopic disordered conduc- 
tors. The 3D conductors are often "neglected" in favor 
of the more popular and tractable playgrounds — two- 
dimensional systems (2D), where one can study states 
resembling 3D critical wave functions in a wide range 
of systems sizes and disorder strengths, 9,10 or quasi one- 
dimensional systems 11 where analytical techniques 12 ' 13 
can handle even non-perturbative phenomena 1 (like the 
ones at small g). For example, in rf = 2 + e di- 
mensions LD transition occurs at weak disorder (weak- 
coupling regime of the corresponding field-theoretical de- 
scription 12 ), while in d > 3 small parameter needed for 



analytical treatment is lacking. In 3D systems critical 
eigenf unctions, exhibiting multifractal scaling, 6 are ex- 
pected only at the mobility edge E c which separates ex- 
tended and localized states inside the energy band. 

The essential physics of disordered conductors is cap- 
tured by studying the quantum dynamics of a non- 
interacting (quasi)particle in a random potential. This 
problem is classically non-integrable, thereby exhibit- 
ing quantum chaos. The concepts unifying disordered 
electron physics with the standard 'clean' (i.e., without 
stochastic disorder) examples of quantum chaos 15 come 
from statistical approaches to the properties of energy 
spectrum and corresponding eigenstates, which cannot 
be computed analytically. While level statistics of disor- 
dered systems have been explored to a great extent, 16 ' 17 
investigation of the statistics of eigenfunctions has been 
initiated only recently. 9 These studies are not only re- 
vealing peculiar spectral properties of random Hamilto- 
nians, but are relevant for the thorough understanding 
of various unusual features of quantum transport in dif- 
fusive metallic conductors (including the ones which are 
proximity coupled to a superconductor 18 ). The standard 
examples are long-time tails in the relaxation of current 19 
or log-normal tails (in d = 2 + e) of the distribution func- 
tion of mesoscopic conductances. 7 Distribution of eigen- 
function amplitudes is found to be relevant for tunnel- 
ing measurements on quantum dots probing the coupling 
to external leads, which depends sensitively on the local 
properties of wave functions. 20 Experiments which are 
the closest to directly delving into the microscopic struc- 
ture of quantum chaotic or disordered wave functions ex- 
ploit the correspondence between the Schrodinger and 
Maxwell equations in microwave cavities. 21 

The study of fluctuations and correlations of eigen- 
function amplitudes 9,22 in mesoscopic systems has led 
to a concept of the so-called "pre-localized" states. 7,19,23 
In 3D delocalized phase, this notion refers to the states 
which have sharp amplitude peaks on the top of an ex- 
tended background. 24 These kind of states appear even 
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FIG. 1. An example of eigenstates in the band center of 
a delocalized phase. The average conductance at half fill- 
ing is g(EF = 0) » 17, entailing anomalous rarity of the 
"pre-localized" states. The disordered conductor is modeled 
by an Anderson model with diagonal disorder Wdd = 4 on 
a simple cubic lattice 12 3 . For plotting of the eigenfunc- 
tion values in 3D, the sites m are mapped onto the lattice 
site numbers £ {1, 1728} in a lexicographic order, i.e., 
m = (m x ,m y ,m z ) h-» 144(m x — 1) + 12(m y — 1) + m z . 



in the diffusive, I <C L < £, metallic (g >• 1) regime, 
but are anomalously rare in such samples (here I is the 
clastic mean free path and £ is the localization length). 
In order to get "experimental" feeling for the structure 
of states with unusually high amplitude spikes, an ex- 
ample is given in Fig. 1; this state is found in a spe- 
cial realization of quenched disorder (out of many ran- 
domly generated impurity configurations) inside the sam- 
ple characterized by large average conductance. Thus, 
pre-localized states are putative precursors of LD tran- 
sition and determine asymptotics of some of the distri- 
bution functions 6 ' 9 studied in open or closed mesoscopic 
systems. In d < 2, where all states are considered to 
be localized, 4 pre-localized states have anomalously short 
localization radius 9 when compared to "ordinary" local- 
ized states in low-dimensional systems. They underlie 23 
the multifractal scaling in weakly localized (g ^> 1) 2D 
conductors of size L smaller than exponentially large 
£ (which plays the role of a phase transition correla- 
tion length 6 £ c = £ in d < 2). In 3D, the correlation 
length £ c (defined as the size of a hypercube for which 6 
g{£,c) = 0(1)) is always microscopic (£ c ~ A^) in good 
metallic samples, and no multifractal scaling is expected. 
The appearance of small regions inside disordered solids 
where eigenstates can have large amplitudes seems to be 
a "strongly pronounced" analog 10 ' 21 of the phenomenon 
of scarring 25 (anomalous enhancement or suppression of 



quantum chaotic wave function intensity on the unsta- 
ble periodic orbits in the corresponding classical system) 
introduced in the guise of generic quantum chaos. 

In general, the study of properties of wave functions 
on a scale smaller than £ should probe quantum effects 
causing evolution of extended into localized states upon 
approaching the LD critical point. In the marginal two- 
dimensional case, the divergent (in the limit L —> oo) 
weak localization (WL) correction 26 to the semiclassical 
Boltzmann conductivity provides an explanation of local- 
ization in terms of the interference between two ampli- 
tudes to return to initial point along the same classical 
path in the opposite directions. 27 This simple quantum 
interference effect leads to a coherent backscattering (i.e., 
suppression of conductivity) in a time-reversal invariant 
systems without spin-orbit interaction. However, in 3D 
systems WL correction is not "strong" enough to provide 
a full microscopic picture of complicated quantum inter- 
ference processes which are responsible for LD transition, 
and facilitate the expansion of quantum intuition. 

The paper presents the statistics of eigenfunction in- 
tensities | \& q, (r) | 2 in isolated 3D mesoscopic conductors 
characterized by two different types of microscopic dis- 
order. Numerical methods employed here make it possi- 
ble to treat phenomena in both semiclassical (described 
by Bloch-Boltzmann formalism and perturbative quan- 
tum corrections 12 ) and fully quantum transport regime 
(dominated by non-perturbative effects, where semiclas- 
sical concepts, like £, loose their meaning 28 ), as well as in 
the crossover realm. Since mesoscopic physics has pro- 
vided efficient techniques 29 for "measuring" exactly the 
transport properties of finite-size samples on the com- 
puter, this study connects the eigenstates statistics of a 
closed sample to its zero-temperature conductance. The 
statistical properties of eigenstates are described by the 
disorder-averaged distribution function 14 ' 23 
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on N discrete points r inside a sample of volume V. Here 
p(E) = (J2 a 6(E — E a )) is the mean level density at 
energy E. Averaging over disorder is denoted by (...). 
Normalization of eigenstates gives t — J dtt f(t) = 1. 
A finite-size disordered sample is modeled by a tight- 
binding Hamiltonian (TBH) with nearest neighbor hop- 
ping integral t mn 
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on a simple cubic lattice 16 x 16 x 16 of lattice constant a. 
Each site m contains a single s-orbital (r|m) = ^(r — m). 
Periodic boundary conditions are chosen in all direc- 
tions. In a random hopping (RH) model the disor- 
der is introduced by taking the off-diagonal matrix el- 
ements to be a uniformly distributed random variable, 30 
1 — 2Wrh < tmn < 1, while diagonal elements are zero 
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FIG. 2. Statistics of wave function intensities in the RH 
Anderson model, with Wrh = 1, on a simple cubic lattice 16 3 . 
The distribution function f(t) in Eq. (1) is computed for the 
states around the following energies: (a) E = 0, (b) E = 1.5, 
(c) E = 2.55, and (d) E = 2.75. Disorder averaging is per- 
formed over JVe„ s = 40 different samples. The Porter- Thomas 
distribution (3) is denoted by RMT. The part of the distribu- 
tions inside the dashed box in enlarged on Fig. 3. 



£m = 0. The strength of the disorder is measured by 
Wrh- The other system studied is described by a diag- 
onally disordered (DD) Anderson model with potential 
energy e m on site m drawn from the uniform distribu- 
tion, — Wdd/2 < e m < VFdd/2, and t mn = 1 is the 
unit of energy. The Hamiltonian (2) is a real symmetric 
matrix because time-reversal symmetry is assumed. The 
results for f(t) in the samples described by the RH and 
DD Anderson models are shown on Figs. 2, 3 and Fig. 4, 
respectively. Although some of the samples are charac- 
terized by similar values of conductance, the eigenstates 
in the two models show different statistical behavior. In 
what follows the meaning of these findings is explained in 
the context of statistical approach to quantum systems 
with non-integrable classical dynamics. In particular, the 
results are compared to the universal predictions of ran- 
dom matrix theory (RMT). 

In the statistical approach 17 of RMT, the Hamiltonian 
of a quantum chaotic system is replaced 31 by a random 
matrix drawn from an ensemble defined by the symme- 
tries under time-reversal and spin-rotation. This leads 
to the Wigner-Dyson (WD) statistics for eigenvalues and 
Porter-Thomas (PT) distribution for eigenfunction inten- 
sities. For the Gaussian orthogonal ensemble (GOE), rel- 
evant for studies of time-reversal-invariant Hamiltonians 
like (2), the PT distribution is given by 
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: exp(-i/2). 



(3) 



The function /pt(^) is plotted as a reference in Figs. 2, 3 
and Fig. 4. The predictions of RMT are universal, de- 
pending only on the symmetry properties of the rele- 



10" 



1(T 



io-' 



10"' 



W =1 

RH 



(A) 

(c) 



L(a)* v, «i 
(bT* 5 ^ coffin 



ili-^RMT 



200 400 600 
t=l¥l 2 V 



800 



FIG. 3. Statistics of wave function intensities in the RH 
Anderson model, with Wrh = 1, on a simple cubic lattice 
with N — 16 3 sites. This Figure plots the same distributions 
f(t) as the ones plotted in Fig. 2, in the range defined by 
the dashed square in Fig. 2. The same labels apply to both 
Figures. 



vant ensemble. They apply to the statistics of real dis- 
ordered systems 33 in the limit g — > oo (g = irE^/A, 
where A = 1/ p(E) is the mean energy level spacing and 
-Etii = KD/L 2 is the Thouless energy, set by the classical 
diffusion across a sample of size L with diffusion constant 
T>). The spectral correlations in RMT are determined 
by logarithmic level repulsion which is independent of 
true dynamics. 17 All sample-specific details are absorbed 
into the mean level spacing 17 A. Also, the level cor- 
relations are independent of the eigenstate correlations. 
The RMT answer (3) for the distribution function (1) 
was derived by Porter and Thomas 34 by assuming that 
the coordinate-representation eigenstate (r|^ a ) in a dis- 
ordered (or classically chaotic system) is a Gaussian ran- 
dom variable. The behavior of ^(r), even within the 
framework of RMT, is simple only in the systems with 
unbroken or completely broken time-reversal symmetry 
[the only difference between the two limiting ensembles 
is the functional form of f{t)}. 35 Thus, RMT implies sta- 
tistical equivalence of eigenstates which equally test the 
random potential all over the sample — typical wave func- 
tion has more or less uniform amplitude 1/yV, up to 
inevitable Gaussian fluctuations. 

Microscopic theory brings corrections to the RMT re- 
sults in the case of samples with finite g. In the finite- 
size systems level statistics follow RMT predictions in 
the ergodic regime, i.e., on the energy separation scale 
smaller than i?Th ■ Non- universal corrections to the spec- 
tral statistics 36 or eigenfunction statistics 9,23,37 (which 
describe the long-range correlations of wave functions) 
depend on dimensionality, shape of the sample, and con- 
ductance g. These deviations from RMT predictions 
grow with increasing disorder (i.e., lowering of g). At 
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FIG. 4. Statistics of wave function intensities in the DD 
Anderson model on a simple cubic lattice 16 3 . The distribu- 
tion function f(t), Eq. (1), is computed for the states around 
following energies. Upper panel, Woo = 10: (a) E — 0, (b) 
E = 6.0, (c) E = 7.45, and (d) E = 7.85. Lower panel, 
Woo = 6: (a) E = 0, (b) E = 4.1, (c) E = 6.56, and (d) 
E — 6.7. Disorder averaging is performed over N^ns = 40 
different samples. The Porter- Thomas distribution (3) is de- 
noted by RMT. 



the LD transition wave functions acquire multifractal 
properties, while the critical level statistics become scale- 
independent. 38 For strong disorder or, at fixed disorder, 
for energies \E\ above the mobility edge \E C \, wave func- 
tions are exponentially localized. The simple (and usu- 
ally invoked) picture is that of a wave function which 
decays as 'F(r) = p(r) exp(— r/£) from its maximum cen- 
tered at some point inside the sample of size L > £. Here 
p(r) is a random function and approximately spherical 
symmetry of decay is assumed. Since two states close 
in energy are localized at different points in space, there 
is almost no overlap between them. Therefore, the lev- 
els become uncorrelated and obey Poisson statistics. If 
p(r) = c is simplified to a normalization constant, the 
distribution function of intensities is given by 
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where a spherically symmetric sample of radius L/2 is 
assumed. In the localized phase £ <C L, f^(t) is expected 
to be insensitive to the assumed shape of the sample, and 
is determined by the ratio of these two relevant length 
scales. 



FIG. 5. Inverse Participation Ratio 1(2) , averaged over 
both 40 different conductors and small energy bins, of eigen- 
states in the RH and DD Anderson models on a simple cubic 
lattice with N = 16 ,i sites. Top: diagonal disorder with (a) 
Wdd = 6, and (b) Woo = 10. Bottom: off-diagonal disorder 
with (a) W RH = 0.25, (b) Wrh = 0.375, and (c) Wrh = 1. 



The distribution function f(t) is equivalently given 
in term of its moments b q — J dtt q f(t) (this state- 
ment is not rigorous since examples of different distri- 
bution functions which possess exactly the same sets 
of moments are encountered in various statistical prob- 
lems 39 ' 40 ). For GOE, the PT distribution (3) has mo- 
ments b^ T = 2W- q+1 T(q + l/2)/r(l/2). They are re- 
lated to the moments I a (q) = J dr \^ a (r)\ 2g of the wave 
function intensity. In the universal regime g — > oo wave 
functions cover the whole volume with only short-range 
correlations (on the scale |ri — r2 1 < i) persisting between 
^ f a (ri) and ^ a (^2)- This means that integration in the 
definition of I a (q) provides self-averaging, and I a (q) does 
not fluctuate in the universal limit, 37 i.e., I a (q) = b^ T . 
On the other hand, at finite g spatial correlations of wave 
function amplitudes at distances comparable to the sys- 
tem size are non- negligible. Therefore, I a (q) fluctuates 
from state to state and from sample to sample. 22,41 Al- 
though these long-range spatial correlations necessitate 
to study the full distribution function 37,41 of I a (q), for 
the subsequent analysis in this study it is enough to use 
an ensemble average of I a (q), i.e., following Wegner 42 



7(g) =A( J] \y a (r)\ 2q 6(E-E a ) 



(5) 



The moment I a (2) is usually called inverse participation 
ratio (IPR) . It is a one- number measure of the degree of 
localization (i.e., it measures the portion of space where 
the amplitude of the wave function differs markedly from 
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FIG. 6. Conductance and density of states in the RH and 
DD Anderson models on a simple cubic lattice 16 3 : RH dis- 
order, (a) and (d), of strength Wrh = 1 (mobility edge is 
at \E C \ — 2.53); diagonal disorder of strength (b) Woo = 6 
{\E C \ ^ 6.5), and (c) Woo = 10 {\E C \ ~ 7.4). Disorder av- 
eraging is performed over A^Ens = 20 different samples for 
conductance and AWs = 40 for DoS. 



zero). This becomes obvious from the scaling properties 
of the average moments I(q) with respect to the system 
size 
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Here d*(q) < d is the fractal dimension. Its dependence 
on q is the hallmark of multifractality of wave functions. 
The multifractal wave functions are delocalized, but ex- 
tremely inhomogeneous occupying only an infinitesimal 
fraction of the sample volume in thermodynamic limit. 
The IPR is affected by mesoscopic fluctuations which 
scale in metallic samples as 22 5I a (2) ~ l/g 2 oc L A ~ 2d . In 
the critical region {g ~ 1) fluctuations 22 are of the same 
order as the average value, which is then not enough to 
characterize the critical eigenstates (even though multi- 
fractal wave function extend throughout the whole sam- 
ple, their IPR is not self-averaging 41 ). 

I use 1(2) as a rough guide in selecting eigenstates 
with different properties in the delocalized phase (Fig. 5). 
The second parameter used in the selection procedure is 
the conductance g(Ep) computed for a band filled up to 
the Fermi energy Ep equal to the state eigenenergy (see 
Fig. 6). The conductance as a function of band filling 
allows one to delineate delocalized from localized phase 
as well as to narrow down the critical region around LD 
transition point (which is defined by g(E c ) ~ 1). Upon 
inspection of these two parameters, a small window is 
placed around chosen energy, and f(t) is computed for 
all eigenenergies whose eigenvalues fall inside the window. 
This provides more detailed information on the structure 
of eigenstates than is encoded in IPR. 

The average IPR for both RH and DD Anderson model 



FIG. 7. Conductance and density of states in the RH An- 
derson models on a simple cubic lattice 16 3 . The strenght of 
the RH disorder is Wnu = 0.25, which is so weak that features 
of a discrete spectrum of a finite-size lattice are still visible 
in the disorder- averaged DOS. Sharp lines correspond to the 
DoS of a clean system on the same lattice (scaled by 1/10 
for clarity). Disorder averaging is performed over Neus = 20 
different samples for conductance and A^Ens = 40 for DoS. 



is shown in Fig. 5. The models with random hopping 
have been around for a long time, 43,44 but have attracted 
considerable attention only recently, inasmuch as they 
show a disorder induced quantum critical point in less 
than three dimensions, 30,45-47 where derealization oc- 
curs in the band center. Furthermore, several models 
used for strongly correlated electron systems share simi- 
lar mathematical structure with the random hopping 2D 
models formulated in the framework of (non-interacting) 
disorder electron physics. 48 Real solids which could be 
described by TBH (2) with off-diagonal disorder include 
doped semiconductors, 44 such as P-doped Si, where hop- 
ping integrals i mn vary exponentially with the distances 
between the orbitals they connect, while diagonal on-site 
energies e m are nearly constant. The behavior of low- 
dimensional RH Anderson model goes against the stan- 
dard mantra of the scaling theory of localization for non- 
interacting systems 4 that all electron states are localized 
in d < 2 for arbitrarily weak disorder. The possibil- 
ity of derealization transition in one dimension 49 goes 
back to the work of Dyson 50 on glasses. Also, the scaling 
theory for quantum wires with off-diagonal disorder re- 
quires two parameters 51 which depend on the microscopic 
model, thus breaking the celebrated universality in dis- 
ordered electron problems. In 3D case explored here, the 
states in the band center are less extended than other 
delocalized states inside the band (Fig. 5). However, the 
off-diagonal disorder is not strong enough 52 to localize 
all states in the band, in contrast to the usual case of 
diagonal disorder where whole band becomes localized 53 
for Wg D > 16.5. 

The mobility edge for the strongest RH disorder 
Hrh = 1, as well as for DD models, is found by looking 
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at an exact zero-temperature static conductance. This 
quantity (which is a Fermi surface property) is computed 
from the Landauer-type formula 54 

g(E F ) = Tr[t(E F )tf{E F )], (7) 

where transmission matrix t{E F ) is expressed in terms 
of the real-space (lattice) Green functions 29,55 for the 
sample attached to two clean semi-infinite leads. To 
study the conductance in the whole band of the DD 
model, f mn — 1.5 is used 55 for the hopping integral 
in the leads. This mesoscopic computational technique 
"opens" the sample, thereby smearing the discrete lev- 
els of an initially isolated system. Thus, the spectrum 
of sample+leads— infinite system becomes continues and 
conductance can be calculated at any E F inside the band. 
However, the computed conductance, for not too small 
disorder 55,56 or coupling to the leads 57 (which are of the 
same transverse width as the sample 56 ) is practically 
equal to the "intrinsic" conductance g = ttEt^/A de- 
termined by the spectral properties of a closed sample. 
The conductance and density of states (DoS) 

N(E) = 2^1, (8) 

are plotted in Fig. 6 for the samples whose eigenstates 
are investigated (the factor of two is for spin degener- 
acy). The DoS is obtained from the histogram of the 
number of eigenvalues which fall into equally spaced en- 
ergy bins along the band. The conductance and DoS of 
the RH model have a peak at E = 0, which becomes 
a logarithmic singularity in the limit of infinite system 
size. 50 To get an insight into the "general weakness" of 
the off-diagonal disorder, Fig. 7 plots g{E F ) and N(E) 
for the low Wrh = 0.25. In this case N(E) still resem- 
bles the DoS of a clean system, even after ensemble aver- 
aging. On the other hand, the conductance is a smooth 
function of energy since discrete levels of an isolated sam- 
ple are broadened by the coupling to leads. The same 
is true for the DoS computed from the imaginary part 
of the Green function for an open system. The mobil- 
ity edge is absent at low RH disorder (Wrh = 0.25 and 
Wrh = 0.375) for system sizes L < 16a. This means that 
localization length £ is greater than 16a for all energies 
inside the band of these systems. For other samples on 
Fig. 6 the mobility edge appears inside the band. This 
is clearly shown for Wrh = 1 case where band edge Eb 
(N(Eb) — 0) differs from E c . The mobility edge is located 
at the minimum energy \E C \ for which g(E c ) is still dif- 
ferent from zero. The conductance of finite samples is al- 
ways finite, although exponentially small at \E F \ > \E C \. 
The approximate values of \E C \ listed in Fig. 6 are such 
that conductance satisfies: g(E F ) < 0.1, for \E F \ > \E C \; 
typically g(E c ) S (0.2,0.5) is obtained, like in the re- 
cent detailed studies 58 of conductance properties at E c . 
Thus found E c is virtually equal to the true mobility 
edge, which is properly defined only in thermodynamic 



limit (and usually obtained from some numerical finite- 
size scaling procedure 2 ). Namely, the position of mobility 
edge extracted in this way will not change 59 when going 
to larger system sizes if £ < L for all energies \E\ > \E C \. 

The distribution /(f) of cigenfunction intensities has 
been studied analytically for diffusive conductors close 
to the universal RMT limit (where conductance is large 
and localization effects are small) in Refs. 9,23 using the 
supermatrix er-model (NLSM), 12 or by means of a "di- 
rect optimal fluctuations method" of Ref. 60. Numer- 
ical studies of statistics of eigenstates in DD Anderson 
model were conducted in 2D for all disorder strengths, 10 
while in 3D the focus has been on the states appearing 
in the semiclassical transport regime where comparison 
with analytical predictions (parameterized by semiclas- 
sical quantities, like k F tj can be made in the regions of 
small 61 and large 62 deviations of /(f) from PT distribu- 
tion. Here I show how /(f) evolves with the strength of 
(different types of) disorder in 3D samples, where gen- 
uine LD transition occurs in the strong coupling regime 
of the corresponding field-theoretical formulation. 

In the weakly disordered (k F £ >• 1) metallic (j > 1) 
conductors pre-localized states are extremely rare. For 
example, in an ensemble of 20 000 samples (Wdd = 4), 
whose typical transport properties are well-described by 
semiclassical theories, only four states would show up in 
the band center which exhibit similar amplitude splashes 
like the one in Fig. 1. Thus, to get a far tail (where de- 
viations from PT distribution are large) of /(f) in such 
conductors one has to search through enormous ensemble 
and locate special configurations of a random potential. 62 
The maximum wave function amplitudes which can be 
observed in this pursuit are, plausibly, determined by the 
strength of disorder, i.e., conductance g. It is, however, 
interesting that tails at small but finite g in the strong 
disorder (like Wdd = 10) can be longer than the tails 
of states in the localized phase, where g is vanishingly 
small, for \E\ > \E C \ at some smaller fixed disorder (like 
in the case of Wdd = 6; compare the two panels in Fig. 4 
using respective conductances from Fig. 6). For strong 
enough disorder long tails of /(f) are found, even by in- 
vestigating small ensembles of disordered conductors (as 
shown below), since the frequency of appearance of pre- 
localized states is greatly enhanced (while system is still 
on the delocalized side of LD transition). 

The complete eigenproblem of a single particle random 
Hamiltonian is solved exactly by numerical diagonaliza- 
tion. Then, /(f) is computed as a histogram of intensities 
for the chosen eigenstates in: delocalized (\E\ < \E C \), lo- 
calized phase (|-E| > |-E C |)) and critical region around the 
mobility edge \E C \. The two delta functions in Eq. (1) 
are approximated by a box function d(x). The width of 
5{E — E a ) is small enough at a specific energy that p(E) 
is constant inside that interval. For each sample, 5-10 
states are picked by the energy bin, which effectively pro- 
vides additional averaging over the disorder (according 
to ergodicity 17 in RMT). The amplitudes of wave func- 
tions are sorted in the bins defined by 5(t — \^ a (r)\ 2 V) 
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whose width is constant on a logarithmic scale. The func- 
tion /(f) is computed at all points inside the sample, i.e., 
N = 16 3 in Eq. (1). 

The evolution of f(t), when sweeping the band through 
the interesting states, is plotted in Figs. 2, 3 for the 
RH disordered sample. Since pre-localized states gen- 
erate slow decay of /(f) at high wave function intensities 
(where PT distribution is negligible), 12 this region is en- 
larged on Fig. 3. This is obvious from the pre-localized 
example in Fig. 1 where state with large amplitude spikes, 
highly unlikely in the framework of RMT, was found in a 
very good metal. The same is trivially true for the local- 
ized states which determine extremely long tails of (f ) . 
Thus, the long asymptotic tails of f(t), appreciably devi- 
ating from PT distribution, are signaling the onset of lo- 
calization. It is interesting that states in the band center 
of RH model, which define the largest zero-temperature 
conductance [g(E F = 0) w 10.2, Var g(E F = 0) w 0.63], 
arc in fact mostly pre-localized. Moreover, both the fre- 
quency of their appearance and high amplitude splashes 
resemble the situation at criticality (NLSM-type calcu- 
lation 47 shows that 3D wave functions, sufficiently close 
to the band center, are always extended for any disor- 
der strength). It might be conjectured that these pre- 
localized states would generate multifractal scaling of 
IPR in the band center. This result, together with the 
DoS and conductance from Figs. 6, 7, shows that phe- 
nomena in the band center of 3D conductors with off- 
diagonal disorder are as intriguing as their much studied 
counterparts in low-dimensional systems. 30 ' 45 ' 46 The ori- 
gin of these phenomena can be traced back to a special 
sublattice, or "crural" 44,63 symmetry obeyed by TBH (2) 
with random hopping and constant on-site energy (lead- 
ing to an eigenspectrum which for E a contains —E a , with 
a special role played by E = 0). In the Wrh = 0.25 and 
Wrh = 0.375 cases all states are extended. Their f(t) 
looks similar to the distribution function for the delo- 
calized states at E = 1.5 in the sample characterized 
by Wrh = 1- The distribution function fe(t) defined 
in Eq. (4) fits reasonably well the numerical f(t) gen- 
erated by the states around E = 2.75, where £ ~ 1.2a 
is extracted for the localization length. Thus, one can 
measure approximately £ in this way 10 even though the 
structure of localized eigenstates can be more compli- 
cated 2 than the simple radially symmetric exponential 
decay used to derive f%(t) [e.g., in the case of DD dis- 
order £ ~ 1.3a is obtained for examined localized states 
in Fig. 4 for Wdd = 10 ensemble, while corresponding 
states in Wdd = 6 ensemble seem to be too close to the 
mobility edge to follow /{(f)]. 

The same statistical analysis is performed for the eigen- 
states of DD Anderson model — a "standard model" in 
the localization theory. Figure 4 plots f(t) at specific 
energies in samples characterized by different conduc- 
tances g(E F = Ei). The conductance g(E = 0) of TBH 
with Wdd = 6 is numerically close to the conductance 
of RH disordered samples with Wrh = 1- Nevertheless, 
comparison of the corresponding distribution functions 



reveals disorder-dependent features 9 which are beyond 
universal corrections (i.e., independent of the details of 
random potential) accounted by the properties of a classi- 
cal diffusion operator 14 ' 37 (the spectrum of — 2?V 2 , with 
appropriate boundary conditions, depends on g, shape 
of the sample and dimensionality; note that conductors 
at W ~ 6 with half-filled band lie on the boundary of 
applicability of such scmiclassical concepts 28 ). To get 
the far tail of the eigenstate statistics in the band cen- 
ter a much larger statistics is needed than used here. 62 
This then makes the observed /(f) in a special case of 
the band center of Wrh = 1 disordered conductor even 
more spectacular because of very large amplitudes found 
in the small ensemble of conductors. Thus, the strong de- 
pendence on the microscopic details of random potential 
demonstrated in this study is somewhat different from 
the disorder-specific short length scale ("ballistic" 9 ) con- 
tributions 9,62 to the standard picture of diffusive NLSM. 
Namely, here it seems that special features of off-diagonal 
disorder in the band center generate completely different 
functional form of the far tail, and not just some disorder- 
specific values of the parameters 9 ' 62 in the exp-log-cube 
asymptotics. 9,23 ' 60 . 

In both models, all computed /(f) intersect PT dis- 
tribution (from below) around 6 < f < 10, and then 
develop tails far above PT values. The length of the 
tails is defined by the largest amplitude exhibited in the 
pre-localized state (like that in Fig. 1). For strong DD, 
Wdd = 10, the conductance g(E F ) is smaller than 3.5. In 
this regime transport becomes "intrinsically diffusive" , 28 
but one can still extract resistivity from the approximate 
Ohmic scaling of disorder-averaged resistance 28 (for those 
fillings where 64 g(E F ) > 2). However, the close proxim- 
ity to the critical region g <~ 1 induces long tails of /(f) 
at all energies throughout the band — a sign of increased 
frequency of appearance of highly inhomogencous states. 
This provides an insight into the microscopic structure of 
eigenstates which carry the current in a non-semiclassical 
transport regime 12 (characterized by the lack of simple 
intuitive concepts, like mean free path I, since unwar- 
ranted use of the Boltzmann theory would give 28 I < a 
in this transport regime although the sample is still far 
away from the LD transition) . 

In conclusion, the statistics of eigenstates in 3D meso- 
scopic conductors, modeled by the tight-binding Hamil- 
tonian on the finite-size simple cubic lattice, have been 
studied. The disorder is introduced either in the poten- 
tial energy (diagonal) or in the hopping integrals (off- 
diagonal). Also calculated are the average inverse par- 
ticipation ratio of cigenfunctions and the exact zero- 
temperature conductance as a function of Fermi energy. 
This comprehensive set of parameters makes it possible 
to compare the eigenstates in nanoscale samples with dif- 
ferent types of disorder, but characterized by similar val- 
ues of conductance. Disorder-specific details, which are 
not parameterized by the conductance alone, are found. 
This is in spite of the fact that dimensionality, shape 
of the sample, and conductance are expected to deter- 
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mine the finite-size corrections to the universal predic- 
tions of random matrix theory, at least in the samples 
which are in the scmiclassical transport regime. The ap- 
pearance of states with large amplitude spikes on the 
of top of RMT like background is clearly demonstrated 
even in good metals. At criticality, the proliferation of 
such "pre-localized" states is directly related to the ex- 
tensively studied multifractal scaling of IPR. However, 
even in the delocalized phase with good metallic prop- 
erties (g S> 1), where the correlation length 6 £ c <~ Xp 
defined by the sample conductance g(£ c ) ~ 1 is micro- 
scopic (L < £ c would naturally account for the multi- 
fractal scaling, 6 like in 2D), pre-localized state with un- 
expectedly high amplitudes are found in the band center 
of random hopping disordered systems. They are inho- 
mogeneous enough to generate extremely long tails of the 
distribution of eigenfunction amplitudes, akin to the ones 
observed at criticality. 
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